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(N| ■ Abstract 

o 

We apply, to the new OPAL data on V+A strange spectral functions, a suitable combination 
of QCD spectral sum rules directly sensitive to the combination of strange quark (m s ) and 
tachyonic gluon (A 2 ) masses . Using the mean value of A 2 from a sum rule analysis of different 
Ph! channels, we deduce the invariant strange quark mass m s = (106tp) MeV to order a^, which 
leads to the running mass m s (2 GeV)= (931 2 !) MeV. We also obtain an interesting though less 
accurate estimate of the CKM angle: \V US \ = 0.217 ± 0.019. 
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1 Introduction 



The determination of the strange quark running mass is of prime importance for low-energy 
phenomenology, for CP-violation and for SUSY-GUT or some other model-buildings. Since 
the advent of QCD, where a precise meaning for the definition of the running quark masses 
within the MS'-scheme P has been provided, a large number of efforts have been devoted to 
the determinations of the strange quark mass 1 using QCD spectral sum rules 2 a la SVZ [5], 
in the pseudoscalar 0, the scalar j7j, the e + e~ jHUH] channels, tau-decay data [TUJ [TTJ [T^] and 
lattice simulations ^3], while some bounds have been also derived from the positivity of the 
spectral functions and from the extraction of the quark condensate [T4"l 13*]. 
In the following, we propose a new method for determining the combination of the strange 
quark mass and tachyonic gluon mass (A 2 ), which will allow us to test the effect of this new 
A 2 -term not present in the standard OPE. This A 2 -term has been introduced in in order 
to mimic the UV renormalon contribution in the resummation of the PT series, where it is 
expected to replace the uncalculated infinite number of terms of the PT series. Its value has 
been estimated from e + e~ into hadrons data [TJJITHj to be: (a s /7r)A 2 ~ —(0.06 ± 0.03) GeV 2 , 
while the pseudoscalar channel and a fit of the lattice data in the x-space for the (pseudo)scalar 
and V+A channels leads to (a s /7r)A 2 ~ — (0.12±0.06) GeV 2 . In the following, we shall consider 
the average of the previous estimates: 

d 2 = (a s /n)X 2 ~ -(0.07 ± 0.03) GeV 2 , (1) 

which is almost scale independent. One should notice that, the effect of this term which is 
relatively negligible does not perturb and, in some cases, improves the well-established existing 
sum rules results [THJH3HH] an d the precise determination of a s from r-decays [171120] obtained 
originally in [21] . It also solves [T7j the sum rule puzzle scales noticed by [22] m the pseudoscalar 
pion and gluonia channels, where the scales are much larger than the one of the p meson. We 
shall further study the effect of d 2 in the determination of the strange quark mass from tau 
decays. In previous analysis of this channel JOIEJE], insphed from the first analysis of tau- 
decay data [21], the A 2 -effect, which is flavour independent, is absent, to leading order of PT, 
as the authors work with the difference of the spectral functions in the ud and us channels (so- 
called flavour breaking sum rules). However, the price to pay is the large cancellation between 
the two different spectral functions, implying a large error bar in the final result. Some other 
eventual problems appearing in these analysis, are the decrease of the value of the strange 
quark mass output (absence of stability in the number of moments) and the deterioration 
of the convergence of the OPE for increasing dimension of the moments, and the non-trivial 
separation of the spin zero and one parts of the spectral functions in some of the analysis. In 
this paper, we avoid these problems by working with a given flavoured us spectral function 
involving the sum of the spin zero and one mesons, but in the same time, our analysis will 
be affected by A 2 . This spectral function has been measured by ALEPH, OPAL and CLEO 
j2H],and more recently by OPAL |24j . 

1 For reviews, see e.g. 
2 For a review see e.g. 0]. 
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2 The sum rules and the QCD expression 



We shall be concerned with the two-point correlator: 

n£ + V) = i Jd 4 x e^(0\TJ,(x) 4(0)|0> = (q»q„ - q 2 g„) H^g 2 ) + gV4°W) (2) 
built from the charged local weak current: 

J» = «7^(! - 7s)s • (3) 

The upper indices 0, 1 refer to the corresponding spin of hadrons entering into the spectral 
function. Following SVZ [S], the correlator can be approximated by: 



d>0 



where 2 d = C 2 d(02d) is the short hand-notation of the QCD non-perturbative condensates 
(0 2 d) of dimension D = 2d and its associated perturbative Wilson coefficient C 2 d', q 2 = —{Q 2 > 
0) is the momentum transfer. The spectral function (v + a): 

-Imll v+A = (v + a) . (5) 
has been measured using r-decay data, via: 

2 ( 2t\- 1 B(r^(V/AY s =-^ J =^ + u T ) 1 dN , 
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M r 2 y 5 (r -> e-p e + i/ T ) N v/A dt 



^B(r^(V/AY s =- 1 ' J =^ + u T ) 1 d tf v/i4 



V o/ a O = r\ir 19 o 1 ~ 77? — 57 7~ 7 \ ~T T , (6) 

±? (r — > e z/ e + v T ) Ny/ A dt 

where OPAL has used \V US \ = 0.2196 ± 0.0023 [25, , M T = 1776.9^;^ MeV gg] and S ew = 
1.0194 ± 0.0040 [23. 

For a pedagogical purpose, we write the QCD expression of the V+A correlator to leading 
order in a s and m s (see e.g. (2D El), and including the leading new A 2 -tachyonic gluon term 

Eh 



with obvious notations, where: 



0\ = 3m 2 + ^A 2 

7T 



O 4 = 47r 2 (m s (^)+m u (^)) + |(a s (G^) 2 ) 
6 ° = 2 -^f* s (uu) 2 (8) 
Its inverse Laplace transform sum rule (LSR) reads 12H1 12U Hj : 

£0 = f <tt e-^-Imn^ = £(l - e-*« - O 2 r + 4 V + :W) , (9) 

JO 7T Z7T Z Z 



from which, one can derive, to leading order in a s and m s , the LSR: 
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& = --^£o = J q dt te-^-Imll^ = — |l - (1 + t c r)e-^ - 0\r 2 - O 6 r 3 | , (10) 
and the leading order FESR [Uli]: 

rU - i- ^(o + D _ u r on 



Mq = 1° dt -ImnStJ 

JO 7T 



27T 2 L 



Ml s / " (iii Ii m n<?S = |j{i-2f]. (ii) 







From these previous sum rules, one can derive the combination of sum rules 

*■ .... . , 



3. 



Mo = C -Td= f C dt (1 - tr) e- tT -ImU^ 

J0 71 



<Sio = M - r-^i = / <c dt f 1 - 2 r) - Imn v°3 > ( 12 ) 

£ c JO V t c / 7T 

which are sensitive, to leading order, to m 2 s and A 2 . Unlike the individual sum rules, these 
combinations of sum rules are less sensitive to the high-energy tail of the spectral functions 
(effect of the t c threshold). Then, one may expect that they are more accurate than the former. 
However, this accuracy may not be comparable with the one of the r-decay like-sum rules where 
threshold effect suppresses completely the effects near the real axis. 
We shall also work with the ratio of moments: 

^io = 2^i, (13) 

which will also be useful for testing the duality between the LHS (experiment) and RHS (QCD 
theory). 



3 QCD corrections and RGI parameters 

In order to account for the radiative corrections, one introduces the expressions of the running 
coupling and masses. 

• To three-loop accuracy, the running coupling can be parametrized as |2T) 
<•.(") = 4 0, {l-«« 0| f loglog^ 

+ ( a <»») 2 | | 1 o g Mo g ^-| 1 o gl o g ^-| + |] + (^)}, (14) 

with: 

f S -Pi log (^/A) (15) 

3 A sum rule similar to A/io has been used for the first time in the pseudoscalar channel for testing the size 
of the SU(3) breakings in the kaon PC AC relation \61\ . 
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and (3i are the 0(a\) coefficients of the /3-function in the MS'-scheme, which read for three 
flavours 0]: 

A = -9/2, & = -8, P 3 = -20.1198. (16) 
• The expression of the running quark mass in terms of the invariant mass rhi is PjJ H] : 



m^v) = m, (-fta„(i/))-» /A |l + j| h- - ^ 0.(1/) 
+ 2 /3 2 U A> i (31 U ft J + A U ft J 



a^(i/) + 1.95168a*(z/) L (17) 



where 7, are the C?(a*) coefficients of the quark- mass anomalous dimension, which read for 
three flavours [I]: 

7l = 2, 72 = 91/12, 73 = 24.8404. (18) 

• To order oP s , the perturbative expression of the correlator reads, in terms of the running 
coupling evaluated at Q 2 = v 2 (see e.g. |2*T| Hj): 



Q ^u v+A m = jl 



1 + K 



a s (Q 2 



+ 1.6398a: + 6.3711o; + .. . 



7T 



(19) 



• The D = 2 contribution reads to order a^, in terms of the running mass and by including the 
A 2 term EH El 11 C3 : 



1 f 



Q 2 U { y-2 ] (Q 2 ) ~ -^2"| a ^ 2 + 3m2 s i 1 + 2-333a s + 19.58a 2 + 202.309^ + Q(a 



(20) 



The coefficient of the a 4 term has been estimated using PT optimisation schemes arguments 
to be K ~ 2276 ± 200. In our approach its effect like all unknown higher order terms will be 
mimiced by the A 2 -term present in the D = 4 contribution given below. 
• The D = 4 contributions read H] : 



q a ^v:1\q 2 ) 



Ylix 



1 - ^a s ) (a s G 2 ) 



+ (1 — a s ) (m u uu + m s ss) + —a s (m s uu + m u ss) 



+ 



4tt 2 



12 -il l 1 — 

— a„ + 1 m 



287T 2 



I 16C(3) ) „„ 



+ ^M-l +4C(3) ) m -°-4 



(21) 



where the last term is due to the A 2 term [T7], and £(3) = 1.202... We shall use as input 
A 3 = (375 ± 25) MeV for three flavours, the value of a s A 2 in Eq. © and ^2 HI I 

(m u + m d ) (uu + dd) = -2m 2 f 'I 
(m s + m u )(ss + uu) ~ — 2 x ^.1m 1 K f 2 K 

(a s G 2 ) ~ (0.07 + 0.01) GeV 4 



pa s (uu)' 



(5.8 ± 0.9) x 10" 4 GeV 6 



(22) 
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where: f n = 93.3 MeV, fx = ^--^f-n and we have taken into account a possible violation of kaon 
PCAC as suggested by the QSSR analysis [HUE]; p — 2 — 3 is the measures of the deviation 
from the vacuum saturation estimate of the four-quark condesates. 

• FESR including radiative corrections can be deduced from the previous expressions of the 
correlator using either the Laplace or the Gaussian transforms |3*TH 0j. For the FESR, the 
most dominant contributions induced by the radiative corrections to the unit operator is the 
dimension-two terms: 

so * = - (£) (I) < • < 23 > 

which adds to the contribution of the terms in Eq. (|20p. One can check that this term remains 
a correction of the A 2 contribution in this previous equation. Some other corrections induced by 
the log- dependence of the running terms taken into account in the numerical analysis remains 
also tiny corrections. 

4 Phenomenological analysis 

We parametrize the spectral function by using the most recent OPAL data discussed in a 
previous section for t until M 2 . In so doing, we parametrize the data using standard polynomials 
fits which delimit the domain spanned by the error bars of the data. At the level of accuracy of 
about 30%, where m s will be determined, the inclusion of the correlations of each data points 
are not necessary. In fact, we have done similar methods in parametrizing the r and e + e~ data 
in the most accurate determination of the g — 2 of the muon, where the error in the estimate 
is at the level of 1% j33]. In this accurate example, our result and the errors agreed quite well 
with ones where the correlations among different data points have been taken into account (33] ■ 
For t above M 2 , we add the QCD step function: 

-lmU v+A (t > M 2 ) ~ 6(t - M 2 )-^J 1 + a s (t) + 1.6398a 2 (t) - 10.284a*(t) 1 , (24) 
7r 2ir z [ J 

which is consistent with the data at t — M 2 . 
Test of duality 

In principle, the value of the t c -cut of the FESR integrals is a free parameter. We fix its 
optimal value by looking for the region where the phenomenological and QCD sides of the ratio 
of moments in Eq. (|T3~|) are equal. We present this analysis in Fig. 1, by showing the value of t c 
predicted by the sum rule versus t c and by comparing the result with the exact solution t c = t c 
for all values of t c . From Fig. 1, one can see that the upper values of the data points provide 
stronger constraints on t c than their central value. Considering this stronger constraint, we 
deduce, that QCD duality is best obtained at: 

i c — M 2 , (25) 
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Figure 1: FESR prediction of t c versus t c . The green curve corresponds the central value of the data; the 
red one to the larger values of the data; the black continuous line is the exact solution t c = t c . The curves 
correspond to the value rh s = (56 — 145) MeV and giving a s X 2 = —0.07 GeV 2 . 



where one expects to get the optimal value of m s from FESR. In order to get this number, we 
have used the value of the invariant mass m s = (56 — 145) MeV, which is a tiny correction in 
this duality test analysis. Once we have fixed the value of t c where the best duality from the 
two sides of FESR has been obtained, we can now estimate some other observables. 

Estimate of rh s versus A 2 

We can, in principle, estimate m s using the combinations of LSR or/and FESR in Eq. (jl2)) . A 
sample of analysis is given in Fig. El for the LSR and in Fig. El from the FESR: 
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Figure 2: LSR prediction of m s in GeV versus the LSR variable r in GeV 2 for -a s \ 2 = 0.07 GeV 2 , using 
the central value of the data and fixing t c = 3.15 GeV 2 . 
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Figure 3: FESR prediction of rh s in GeV versus t c in GeV 2 for -a s A 2 = 0.07 GcV 2 . The green curve 
corresponds to the central value of the data; the red one to the larger values of the data. 



Table 1: Estimated value of rh s versus A 2 



-a s X 2 in GeV 2 rh s in MeV 



0.02 


53 ± 38 ± 8 


0.04 


79 ± 29 ± 8 


0.06 


98 ± 24 ± 9 


0.07 


106 ± 23 ± 9 


0.08 


114 ±21 ±9 


0.10 


128 ±19 ±11 


0.12 


140 ±18 ±11 


0.15 


157 ±16 ±11 



The first error is due to the data, the second one to A. The error due to the choice of t c around 
the duality region is quite small of about 3 MeV as can be seen in Fig. |21 



• In Fig. we give the value of rh s versus the LSR variable r at given t c = M 2 . One can 
notice that, though the "optimal" estimate looks reasonnable, it is obtained at the value 
of r of about 1 GeV -2 , where, at this scale, the PT series of the m 2 coefficient behaves 
badly rendering this result inaccurate. Therefore, we do not consider this result in our 
final estimate. One should notice that the size of this optimization scale is typical for the 
LSR, as the exponential amplifies the low energy contribution to the spectral function. 
This is not the case of the FESR which acts in the opposite region of the spectral function. 



• In Fig. |2J we give the value of rh s from FESR versus the t c variable. Optimal estimate 
is obtained for the t c - value where the duality between the two sides of FESR is the best. 
One can notice that unlike the case of LSR, FESR results are obtained at t c = 3.15 GeV 2 , 
where the PT series in the m 2 make sense. 

• Results for different values of A 2 in the range given in Eq. (0) are shown in Table ^ 
Considering the mean value of A 2 from Eq. (flj). we deduce from Tabled the estimate: 

m s ~ (106ip) MeV for a s A 2 ~ -(0.07 ±0.03) GeV 2 . (26) 

Using the relation in Eq. (|17|): 

m s (2 GeV) ~ 0.876 m s , (27) 
we translate the result on m s into the value of the running mass at 2 GeV to order a^: 

m s (2 GeV) ~ (93±$j) MeV . (28) 

• This result can be compared with the recent determinations from flavour breaking-sum 
rule in tau-decays [HJQ2], which is not affected to leading order by A 2 : 

m s (2 GeV) ~ (81 ± 22) MeV , (29) 

and with the one: 

m s (2 GeV) ~ (105 ± 26) MeV , (30) 
deduced from the pion sum rule to order a]j, where the A 2 term has been added [T7| l3~| 14] : 

(m u + m d )(2 GeV) ~ (8.6 ± 2.1) MeV , (31) 

plus the ChPT ratio (m u + md)/2m s = 24 .4 ±1.5 [2]. From the previous values, we can 
deduce the average 4 : 

(m s {2 GeV)) ~ (92 ± 15) MeV , (32) 
which is comparable with lattice determinations. 

• Finally, one can compare this result with the lower bounds obtained from the positivity 
of spectral functions proposed in [T3j and updated to order to be (71.4 ± 3.7) MeV in 
[3] (the inclusion of A 2 decreases this value by 5%), and from an independent lower bound 
of about 80-90 MeV derived in ^H] from direct extractions of the quark condensate. Our 
result in Eq. l2*8lis compatible with these bounds. However, values of m s corresponding to 
—A 2 < 0.04 GeV 2 are less favoured by these bounds and by the previous estimates from 
other channels. 

4 We have not included in the average the value of m s from the vector current in |3fi| . which is now under 
reconsideration. 
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Estimate of \V US \ 

We can estimate the CKM angle by working with the FESR M.q or/and the LSR Cq. FESR 
gives the result: 

\V ua \ ~ 0.215 ±0.017 , (33) 

while the LSR gives: 

|K, | ~ 0.219 ±0.021 . (34) 

One can notice that the result is not sensitive to the change of m s in the range 56 to 145 MeV, 
which, a posteriori justifies the use of the exponential Laplace sum rule (LSR), though, in the 
optimal region, the PT series of the m 2 s contribution behaves badly. We take, as a final result, 
the arithmetic average of the FESR and LSR results: 

|K.| ^ 0.217 ±0.019 , (35) 

which, despite the large error, is an interesting output per se as it shows the consistency of the 
approach used to get m s . 

5 Conclusions 

Our result for the strange quark mass, in Eq. ([28)1 . from the strange spectral function of hadronic 
tau decays shows that the one for non-zero value of the tachyonic gluon mass is in better 
agreement with the existing determinations of this quantity from r-decays and pseudoscalar 
channels and with the lower bounds derived from different sum rules EE El E] than the one 
in the theory without a tachyonic gluon. 

Our result in Eq. for V us is less accurate than existing determinations but interesting per 
se for testing the consistency of the whole approach. 

Improvments of the results obtained in this paper require more accurate data, which will also 
help for a much better determination of m s or/and of the tachyonic gluon mass A 2 . 
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